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Abstract. We extend the invariant manifold method for analyzing the asymptotics of dissipative partial differential 
equations on unbounded spatial domains to treat equations in which the linear part has order greater than two. One 
important example of this type of equation which we analyze in some detail is the Cahn-Hilliard equation. We analyze 
the marginally stable solutions of this equation in some detail. A second context in which such equations arise is in the 
Ginzburg-Landau equation, or other pattern forming equations, near a codimension-two bifurcation. 
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1. Introduction and statement of results 

In this paper, we extend the methods developed in [Wl], [W2], [EWW], to study the asymptotic 
behavior of marginally stable non-linear PDE's. These are PDE's such as 

d^u = P{-iV,)u + W\u), 

where u = u(x,t), with x e R'', and where P is a polynomial. In the papers cited above, we 
have treated essentially parabolic problems, i.e., the case where -P(0 = this paper, 

we extend the problem to non-parabolic cases such as P{^) = where P{—iV^) has 

continuous spectrum all the way up to 0. We deal in particular with the stability analysis of 
the Cahn-Hilliard equation [CH] in an infinite domain. Where appropriate, we indicate how to 
formulate the assumptions for more general differential operators and non-linearities. 

The Cahn-Hilliard equation models the dynamics of a material with the following 3 prop- 
erties: 

i) The material prefers one of two concentrations that can coexist at a given temperature. 

ii) The material prefers to be spatially uniform. 

iii) The total mass is conserved. 

The first point above means that we should consider a potential with 2 minima with equal 
critical values, and for concreteness, we will choose W{u) = (1 - u^f.* The Cahn-Hilliard 
equation is then 

= A{-/\u + W'{u)) , (1.1) 

or, expanding, 

d^u = -A^u - 4Au + 4Au^ . (1.2) 
We will be interested specifically in the non-linear stability of the spatially uniform states, 

u{x^ t) = Uq. 

It is obvious that constants are solutions of Eq.(1.2), for any Uq. Furthermore, it is easy to 
check that these solutions are (locally) linearly stable for \uq\ > 3~^^^, and linearly unstable for 
\uq\ < 3~^/^. We concentrate our analysis on the remaining case, namely Uq = ±3"^^^. In this 
case, linearizing about Uq = 3~^^^ leads to the linear equation 

d^v = -A\ , (1.3) 

which has spectrum in (— oo, 0] and corresponds to the case P{^) = It is obvious that 
bounded initial data lead to solutions which tend to as t oo and the purpose of this paper is 
to study under which conditions the addition of the nonlinear terms does not change the stability 
of the solutions. This is more difficult, for two reasons: First, as we have said, the spectrum of 
the linearized problem extends all the way to 0, and second, the nonlinearity does not have a 
sign. 



* In our example, the curvatures of the two minima are equal. This does not seem to be necessary for our proofs. 
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Another, more complicated, example of a similar nature is provided by time-independent 
solutions of the Ginzburg-Landau equations (on R) 

d^u = dl.u -\- u — u\u\^ , (1.4) 

which are exactly on the borderline between being Eckhaus stable and Eckhaus unstable. These 
solutions are 

with q = 1/^/3, cf. [EG]. We will not prove that this problem scales like the Cahn-Hilliard 
equations, but only describe a program which we believe would lead to a proof. The first part 
of the analysis of this problem would follow rather closely that given in [EEW] for the Swift- 
Hohenberg equation. Letting u* = for the critical value q = \ / V^, and writing u = u* + v, 
the equation for v is 

^^v = dlv + v- 2v\u*\^ - v{u*f + 0{v^) . (1.5) 

It has a linear part which is like a Schrodinger operator in a periodic potential (the inhomogeneity 
u*). This can be handled by going to Floquet variables, namely setting 

dke^'"'v,.{x,t) , 

-q 

where is Tr/q'-periodic in x: 

v,{x,t) = J^e^^-'^XmW- 

meZ 

The linear part of Eq.(1.5) leaves the subspace spanned by the invariant, and has discrete 
spectrum in each such subspace. The spectrum is in cr < and the largest eigenvalue is —0(k'^) 
when q equals its critical value q = l/\/3 (which is the case we discuss here). In this sense, 
the problem of the marginal Eckhaus instability resembles the problem of the Cahn-Hilliard 
equation. At this point, the discussion of the problem follows the techniques we developed 
in [EWW]. We would like to rescale as we will do below for the Cahn-Hilliard equation and 
its generalizations, but the problem will be more complicated because the Brillouin zone is 
restricted to A; G [— q]. We then have to check that the non-linearity is "irrelevant" in the 
terminology developed below. Again, as in [EWW], we believe that this will not be quite the 
case, but the saving grace will be that the projection of the potentially non-irrelevant modes 
onto the eigenstates corresponding to the — 0(A;'^) term vanish to some higher degrees because 
of translation invariance of the original problem, cf. [EWW, Section 4], and [S]. 

We place our examples in the following more general setting. Consider equations of the 

form 

(/I I 

— = (-l)-+iA-« + F(«,{a»), (1.6) 
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where the muhi-indices a satisfy \q\ < In — 1, and x G R'^, t > \. Furthermore, F is a 
polynomial in u and its derivatives. We wish to study the asymptotics of the solution u of (1.6) 
as t — > oo. First, one introduces scaling variables by defining 

1 X 

^(^'*) = ^d7(2^ ^(^17(2^1) 'l^sO • (1-7) 

Introducing new variables ^ = xjt^l^'^'^^ and r = log(t + tg), with an arbitrary positive 
constant, the Eq.(1.6) is transformed to the non-autonomous problem 

1^ = (-1)"+! A-^ + . V^v +^v + e(^)-F(e-^^, {e-^^)-^"^}) . (1.8) 

The analysis of this equation involves two steps: 

i) An analysis of the linear operator 

ii) A determination of which non-linear terms are relevant. 

As we will see, the term 1 / (2n)^ ■ plays an important role in the analysis of this linear 
operator as it allows us to push the continuous spectrum of the operator more and more into the 
stable region by working in Sobolev spaces with higher and higher polynomial weights. These 
weights force the functions to decrease more and more rapidly near |a;| = oo. Taking Fourier 
transforms on both sides of Eq.(1.8) we obtain: 

= _(p.p)-i)_ p. Vp{) + e^^)^F*(e-^{;,{e-(^)"(-zp)"{)}) , (1.9) 

where F* is the polynomial F, written in terms of convolution products, (see the discussion of 
the non-linearities below). 

We will discuss the form of the non-linear terms below, and consider first the linear operator 

£ = -{p.pY-l-p.v^. (1.10) 

A straightforward calculation shows that L has the countable set of eigenvalues 

A, = -^, j = 0,l,2,..., (1.11) 
with eigenfunctions (written in multi-index notation), 

= p"e-(^-^)" , (1.12) 

and \a\ = j. 

If we consider L as acting on the Sobolev spaces 

He,m = : lb"^p^llL2 < , for all |a| < £, \(3\ < m} , 
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then XL will have continuous spectrum in the half-plane Re A < —cy^ in addition to the eigen- 
values above. By choosing m appropriately, we can force this continuous spectrum arbitrarily 
far into the left half-plane, and the dominant behavior of the linear operator will be dictated by 
the eigenvalues with the largest real part. 

Remark. In order to switch back and forth from the Fourier transform representation of L 
to the un-Fourier transformed representation of this operator with ease, we also consider the 
Sobolev spaces 

= : \\d>'^^\\u < oo , for all \a\ < £, \p\ < m} . 

Note that Fourier transformation is an isomorphism from ^io ^. 

Note that L is not sectorial, and therefore we know of no way to bound the semi-group 
generated by L by spectral information alone. However, in Appendix A, we develop an integral 
representation of the semi-group and we then show that it satisfies the estimates needed for the 
invariant manifold theorem. 

We next discuss which terms in the non-linearity are "relevant." Consider a monomial 

A-i[{dfuf , (1.13) 

where the a^^^ are distinct multi-indices. After rescaling and taking Fourier transforms this 
becomes 

X ((-ip)"^%)*'% ■ ■ ■ * ((-^p)-'^'^)*'^ 

Here, * denotes the convolution product. If we combine the powers of r in the exponential, we 
see that if 

s 

2n + d < J2{\(^^^^\ + d)kj , (1.15) 
i=o 

then the coefficient of this term will go to zero exponentially fast in r, and hence it will be 
irrelevant from the point of view of the long time behavior of the solutions. 

Definitions. A monomial like (1.14) is called irrelevant if it satisfies the inequality (1.15). It 
is called critical if the l.h.s. of Eq.(1.15) is equal to the r.h.s, and relevant in the remaining case. 
These definitions are suggested by the following which is our first main result: 

Theorem 1.1. Assume all terms in the non-linearity in Eq.(1.6) are irrelevant. For any solution 
u{x, t) of Eq.(1.6) with sufficiently small initial conditions in ^ (with £ > {In — 1) + d/1 
and m > {n + 2)/ (2n) ), there is a constant B* , depending on the initial conditions, such that 
for every e > 0, 



lim t^^-'^) 

t—>oo 



B* ^ X 



= 0. 
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Here, 

= J2^^ J ^'P^''-'^-^'-''" ■ (1-16) 

Remark. This theorem is a special case of a more detailed analysis which will be given below. 
That analysis will allow us to compute, in principle, the form of the solutions of Eq.(1.6) up 
to 0(t~^), for any A; > 0. We note that if one only wanted the first order asymptotics of the 
solution, one could also use the renormalization group analysis of [BKL]. 

We now apply the Theorem 1. 1 to the Cahn-Hilliard equation. Writing u = 3~^^'^ + w, the 
function w is seen to satisfy 

^ = -A^w + V3A{w^) + Aiw^) . (1.17) 

Upon expanding A{w'^) we obtain two types of terms — those of the form w{dl.w) and those of 
the form (d^.w)^. In both cases, 

J2i\(^^'^\+d)kj = 2d + 2. 

Since n = 2 in this example, these terms will be irrelevant if4+d<2d+2, that is in dimensions 
d > 2. Also, the term A(w^) is irrelevant for d > 1. Thus, as a corollary to Theorem 1.1 we 
get immediately 

Corollary 1.2. Solutions of the Cahn-Hilliard equation in dimension d > 3, with initial 
conditions sufficiently close (in j) to the constant solution u = 3~^^^ behave asymptotically 
as 

= ^ + (^) + ^iw^^i • ^'-'^^ 

Remark. We will examine below what happens in the cases d — 1,2. The case (i = 2 is of 
particular interest because its non-linearity is critical in the renormalization group terminology. 



2. Invariant manifolds 

Note that spectral subspaces corresponding to eigenvalues of £ are automatically invariant 
manifolds for the semi-flow defined by the linear part of Eq.(1.9). The aim of this section is to 
demonstrate that the full non-linear problem has similar invariant manifolds in a neighborhood 
of the origin. This then shows that the conceptual understanding of what is happening can be 
gained purely from a knowledge of (and the scaling behavior of the non-linearity). 
We begin with a proposition concerning the linear semi-group generated by XL. 

Proposition 2.1. Let denote the projection onto the spectral subspace associated with the 
eigenvalues { ^ }j^o' Qk — (^^ ^e,m)- ^^"^ > {n + k + l)/(2n), then there exists 

> such that the semi-group generated by £ satisfies 

\\Qke^^QkV\km < SkM^^)Mi-,,m,Q = 0,h---,^n-l. (2.1) 



Cahn- Milliard Equation 



7 



Proof. The proof, which is presented in Appendix A, is modeled on the proof in [EWW] which 
treats the case n — I. 

Given such estimates on the Unear evolution, the construction of invariant manifolds is 
straightforward. Denote by y the coordinates on the (finite-dimensional) range of P^., and let 
z = Qi^v. Finally let rj — e~'^l^^'^^ — {t-\- ^0)"^/*^^"-*. Then, applying the projection operators 
and Qa; to Eq.(1.9), it can be written as the system of equations 

y = ^ky + fiy^^^v) , 

z = Qj^Lz^ g{y,z,r]) , (2.2) 

V = -^v, 

where ' denotes differentiation w.r.t. r. We next need a bound on the non-linearity: 
Lemma 2.2. Assume u e £ with £ > 2n — I + d/2, and assume 

s 

2n + d < + d\kj . 

Then the non-linear tenn Eq.(1.14) has ^^r„,^_2n-i-i ^orm bounded by 

s 

n ll"llm,^-2n-|-l = '^^ll'*^llm,^-2n-hl ' 

i=o 

withp= Zj=o + d\kj - {In + d)andK = ^ -=0 

Proof. Taking the inverse Fourier transform of Eq.(1.14), and substituting rj — e^'^^^^'^\ 
Eq.(1.14) becomes 

j=0 

The result then follows from the Sobolev embedding theorem. Note that the lemma has the 
immediate corollary (because F is a polynomial): 

Corollary 2.3. Under the hypotheses of Lemma 2.2, for every r > I, the non-linear term in 
Eq.(1.9) is a function from R x g to e-2n+\ ■ 

This corollary in turn implies that the terms in Eq.(1.14) and Eq.(2.2) are all functions. 
This, in conjunction with the estimates on the linear semi-group is sufficient to establish the 
following 

Theorem 2.4. Suppose thati > In — \ -\- d/2 and m> {n + k + I)/ {2n). Suppose further 
that all terms in the nonlineadty satisfy 

s 

2n-\-d < + d\kj . (2.3) 
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Then there exists a C^"*"" function h{ri, y), with a > 0, defined in some neighborhood of the 
origin in R x R'^'™ range(Pfc) ^ ^^^j^ ^^^^ manifold z — h{r], y) is left invariant by the semi-How 
of Eq.(2.2). Furthermore, any solution ofEq.(2.2) which remains near the origin for all times 
approaches a solution of (2.2) — restricted to the invariant manifold — at a rate (e ] . 

Proof. The existence of the invariant manifold, given the assumptions on the linear semi-group 
and the non-linearity, seems, to our knowledge, not to be explicitly spelled out in the literature. 
The formulation which comes closest to our needs is the one given in [H], where the assumptions 
on the non-linearity are those we have in our case, but the semi-group is supposed to be analytic. 
However, Henry's construction of the invariant manifold only uses certain bounds on the decay 
of the semi-group, and not the stronger assumption of analyticity. Those bounds are true in our 
case, by Lemma 2.2. Thus, existence of the invariant manifold follows in fact from Henry's 
proof. 

Once one knows that the manifold exists, it is also easy to show that any solution near the 
origin must approach a solution on the invariant manifold (see, e.g. [C]). Note that even though 
our non-linearity is quite smooth, we cannot hope, in general, to obtain an invariant manifold 
whose smoothness is greater than C^^", since this smoothness is related to the gap between the 
spectrum of A^, and that of Qj^LQ^^, (see, e.g. [LW]). 

3. Applications 

Here, we show how the existence of the invariant manifold implies Theorem 1.1 and related 
results. To prove Theorem 1.1, we assume that all terms in the non-linearity are irrelevant. 
This means that Eq.(2.3) holds. Suppose further that k = \ and that I > In — \ + d/1 and 
m > {n + 2)/{2n). These hypotheses guarantee that Theorem 2.4 applies and hence any 

solution near the origin must approach a solution on the invariant manifold, at a rate (e 2" j 

in H^^g. 

The equations on the invariant manifold can be written as a system of ordinary differential 
equations: 

vo = {'Po\f{y^Hv,y),v)) , 

= + ('^ul/(2/'M^'2/)'^)> ' i=i'--->c?> (3.1) 

V = -j^v, 

where <^q and (flj are the projections onto the spectral subspace of Ag and Aj = — l/(2n), 
respectively. Note that Aj has a cZ-dimensional spectral subspace. 

The important observation to make at this point is that since the non-linearity is assumed 
to be irrelevant, there exist constants Cq and Cj such that 

\{^;\f{y,h{r),y),ri))\ < C,rj\ \{^l\f{y,h{r),y),r^))\ < C,rf , 

forsomep > 1. Since ?7(t) = e~^/^^"^?7(0), this implies immediately that solutions of Eq.(3.1) 
behave as 
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The eigenfunction with eigenvalue of XL is e or taking inverse Fourier transform, /*, 
cf. Eq.(1.16). Thus, in -ff^ ^ solutions of Eq.(1.9) behave as 

Reverting from scaling variables to the unsealed variables u{x, t) and using the Sobolev lemma 
to estimate the L°° norm in terms of the ^ norm, we obtain Theorem 1.1. Since we 
observed above that the non-linearity in the Cahn-Hilliard equation is irrelevant when d >3, 
we immediately see in this case that Eq.(1.18) holds for initial conditions which are close to 
u = 3 ' , which yields Corollary 1.2. 

4. The critical case 

We now consider the Cahn-Hilliard equation in dimension d = 1, which is the critical case 
in terms of the renormalization group terminology [BKL]. This means that in some non-linear 
terms the inequality Eq.(1.15) becomes an equality. 

In the Cahn-Hilliard equation, when d — 2 (and n = 2), we see that the quadratic term is 
critical, and the cubic term is irrelevant. Note that Theorem 2.4 still implies the existence of an 
invariant manifold tangent at the origin to the eigenspace of Aq. This means that when written 
in the form of Eq.(2.2), the non-linearity can be written as the sum of 2 pieces — one quadratic 
in y and z which is independent of rj (and hence critical) and a cubic piece in y and z which is 
linear in rj (and hence irrelevant). This implies that the Eqs.(3.1), when reduced to the invariant 
manifold, take the form 

% = {vl\f^^\yMv.y).v)) + {^l\f'-'\y.h{Ti,y),rj)) , 

yij = -iyij + <^ul/^'n2/'%'y)'^)> + fel/^'n2/'%'y)'^)>' i = i,2, (4.1) 
V = -iv ■ 

We now exploit the form of the non-linear term in Eq.(l. 17), namely 3^^^A{w^) + A{w^), 
plus the fact that the eigenfunction cpQ = I. Thus if we integrate by parts, we find that 

{^o\f^^Ky,Krj,y),rj)) + {^;\f^'\y,h{ri,y),rj)) =0, 

so that in Eq.(4.1), i/q = and thus j/o(*) = yoi^)- reflects the fact that the Cahn-Hilliard 
equation conserves mass. 

Since from Eq.(4.1) we also see that j/j ^ = 0{e~'^^'^), we find upon reverting to the 
unsealed variables the second main result: 

Theorem 4.1. For d = 2, if the initial conditions of the Cahn-Hilliard equation are sufficiently 
close in H^ ^ to the stationary state u = 3"^^^, then the solution behaves asymptotically as 
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5. The relevant case 

Here, we consider the case of d = 1 where one term of the non-Unearity is relevant. This 
necessitates a change of strategy, because the quadratic term is proportional to r]~^ and hence 
the non-linear terms in Eq.(2.2) are not smooth enough to apply the invariant manifold theorem. 
In order to circumvent this difficulty, we choose a scaling different from Eq.(1.7). Consider 
again the Cahn-Hilliard equation, Eq.(1.17), with u = 3"^^^ + w. Ind = I, we get 

Now let w{x, t) = r^l^Wixjt^l^, logt). Then W satisfies 

^ = -d\w + \i ■ d^w + + 3i/2a|(H^2) + e-^i''d\{yv'') . (5.2) 

Proceeding as in the other cases, we define the linear operator Lj = — + + which in 
Fourier variables becomes 

so that it has eigenvalues [i^ = j = 0, 1, . . . . Thus, unlike the operator L, we have one 
eigenvalue lying in the right half-plane. Let fj = e~'^l^ , and let and denote the amplitudes 
of the eigenvectors with eigenvalues \Xq and /Xj. Then Eq.(5.2) takes the form 

% = i% + /o(%'2/i,^,y^) , 
vx = fi{yo,yi,v,y^) , 
V = -b^ 

y^ = QL^y^ + f^{yQ,yi,r],y^) . 

Here, Q is the projection onto the complement of the eigenspaces corresponding to /Xq and 
y~^ — QW, and /q, /j, and /"""are the projections of the non-linearity onto the various 
subspaces. 

Since the spectrum of QL^Q lies in the half-plane Re // < | , we can construct an invariant 
manifold for Eq.(5.3) which is the graph of a function h^{yQ,y^,ri), and every solution of 
Eq.(5.3) which remains in a neighborhood of the origin will approach this manifold at a rate 
0(e-^/4). What is more, the equations on the invariant manifold are extremely simple in this 
case, since the projections onto the "0" and "1" components correspond to integrating with 
respect to the functions 1 and x, respectively. Applying these projections to the non-linearity 
and integrating once, resp. twice by parts, we see that these projections of the non-linear terms 
vanish. Thus, the equations on the invariant manifold of Eq.(5.3) are simply 



% = 3% > 2/1 = 0, r) = -i?7 
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Thus, as long as the solution of Eq.(5.3) remains in a neighborhood of the origin, it will be of 
the form 

Voir) = e^/\o(0) , 
ffir) = e-^"m 



- .-VS.,.. (5-4) 



Thus, we see that the solution either leaves the neighborhood of the origin, or its asymptotic 
behavior can be read off from Eq.(5.4). Note that the solutions that remain near the origin must 
have Uq = 0. Thus: 

Theorem 5.1. Suppose that the initial condition of the Cahn-Hilliard equation is of the form 
Uq = 3~^/^ + Wq with Wq small in the H^ (^ norm for some £ > 4 and m > 2. Assume 
furthennore that dx Wq{x) = 0. Then the solution is of the form 



where 



dxxw,{x) , r*(e) = ^ / dpe'^^e-^ . 
-co v27r J-oo 

Remark. The constant B* in Theorem 1.1 is not as easy to describe because there, the 
non-linearity in the equation for j/q did not necessarily disappear. 

Proof. The proof is an obvious modification of the one of Theorem 1.1, taking into account 
the special form of the eigenfunctions corresponding to the eigenvalues /Iq and //j. 



Appendix. Bounds on the linear semi-group 

In this appendix, we sketch the proof of Proposition 2.1. The proof is quite similar to the 
estimates on the linear semi-group in Appendix B of [EWW], (which was given for the case of 
a one-dimensional Laplacian, or in the present notation n — d = 1) so we concentrate only on 
the points where the present argument differs from the one in [EWW]. 
We begin with the representation 

Td 

{e^''){x) = 0/ d''zg{z,T)v{e^{x + z)) , (A.l) 

where 

g{z,T) = j d^/!;e^'='^exp(-(A;-A;)"(l-e-")) . (A.2) 

As in [EWW], the action of the semi-group is analyzed by considering separately the behavior of 
the part far from the origin and that close to the origin. The new difficulty here is that we do not 
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have an explicit representation of g as in the case n = 1 . However, the technique of estimating 
the long-time behavior will remain essentially the same. Let Xr be a smooth characteristic 
function which vanishes for |a;| < R and is equal to 1 for | a; | > 4R/3. We start by studying the 
region far from the origin. The analog of Proposition B.2 of [EWW] is 

Proposition A.l. Foreveryi > 1 andeverym > 0, there exist a > and aC{£,m) < oo 
such that for all v e H^ ,^ one has 

(a(T))^" ^ ^ 

for = 0, 1, . . . , 2n — 1. Here, a(T) = 1 — e~^. 

The crucial step in proving this estimate is to derive the asymptotics of g{z, r) for large 
z. This will replace the explicit (Gaussian) estimates for the d = I, n — 1 case analyzed in 
[EWW]. This estimate is provided by the following 

Proposition A. 2. The kernel g{z, r) decays faster than any inverse power of z for \z\ large. 
In fact, one has the estimate 

\giz,T)\ < Ca(T)-^exp(-7(|^|27a(T))5^) , (A.4) 

for some 7 = 7(n, d) > 0. 

Remark. If n = 1, we recover the explicit bound on the Green's function: 



Proof. We need to estimate the quantity 

By rotational symmetry, it suffices to bound the preceding expression for x lying on the positive 
real axis. Setting x = 2na(T)z^'^~^, and k = (p, q), with p e R, and q e R'^"^ this means 
that we must bound 

X = y dpd'^"^gexp(-a(T)(/ + g • g)" + 2inpa{T)z^''~^) . 
If we rescale the variables asp — zt, and q — zs, then we have 

X = z"^ I dtd'^-hex^{-a{T)z^'' {{t^ + s ■ + lint)) . 
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Remark. Note that the polynomial {t^ + s ■ s)'^ + lint is independent of z. 

We will bound X by taking advantage of the fact that the integrand is an entire function 
and translate the contour of integration so that it passes through at least one critical point of the 
exponent. These critical points occur at s = and the roots of t^'^~^ = i - i.e., at the points 
= exp(ifg±J}), = 0, 1 , 2, . . . , 2n - 2. 

Inserting this expression into the exponent of the integrand of X, we see that the value of 
the polynomial at the critical points is 

-a{r)z (exp(^^^^--^)-2mexp(^^^^--^) 

,N / N 2n ,.7r(4/c + 2n)^ 
= (2„-l)a(r).^"exp(.^p^). 

In particular, if we take k = 0, then the real part of the critical value is 

(2n - l)a(T)^2" cos(7r/2 + 7r/(4n - 2)) ^ -a{T)z^'' ■ ^ , 

when n is large (and is negative for all n > 0). Integrating over the region R + observing 
that there is only one critical point on this line, we get, using the techniques of Hormander: 

X ^ a(T)'^/(2")e-^"«(^)^" , 

with C > and C„ — > 7r/2 as n ^ oo, when z ^ oo. Reverting to the original variables, this 
leads to 

where = C^/(2n)^/(2n-i)_ 

We now consider the action of the semi-group on functions localized inside a ball of radius 
R. A key observation here is the following lemma. Let v'o(^) denote the eigenfunction of L 
with eigenvalue and let T{x) = 0q{xY^^ (note that 0q{x) > for all x). 

Lemma A. 3. The operator H — T~^LT is self-adjoint on (a dense domain in) L^(R^) and 
has the same eigenvalues as L. 

Proof. The proof is a straightforward calculation. We note further that if (p^ are the eigen- 
functions of L then the eigenfunctions of H are ipQ,{x) = (T~ Va) (^)- 

If we take the inverse Fourier transform of the eigenfunctions (f^{p) of Eq.(1.12), we see 

that 

\(fi^{x)\ C|a;|'°'' exp(— 7|a;|2^) , 

for some 7 > 0, using the same sort of estimates as those used to bound the kernel g of the 
semi-group. Thus, for |a;| sufficiently large, we get 

\i'ai^)\ ~ exp(— i7|a;|2^) . 
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The usefulness of introducing the operator H is that it is sectorial, since it is self-adjoint 
and bounded below. Therefore, the associated semi-group can be estimated from spectral 
information alone. In particular, if P^. denotes the projection onto the spectral subspace spanned 
by the eigenfunctions with eigenvalues 0> — ik' In '• • • ' In ' Qk^^ defined by = I — Pj^, 
then we have a bound on the operator norm of Qj^e^^Q^. 

HQke^'^Qkll < C,e--^-/(2n) _ (A.6) 

We can use this information to bound the semi-group associated with L. Note that if we denote 
by P^"^ and Q^^ the projection associated with the spectral subspaces of L (as we did for H), 
the we have the identity: 

e^-Ql% = e^-Qf (1 - X> + e^-Qf X^. . 

Since x^f is localized away from the origin, it can be studied with the help of Proposition A.l, 
so we only focus on the other term. And there we get 

lle^^gf (l-Xfi)^ll,,. = \\TT-'e-'^TT-'Q'~^^TT-'{l-XR)v\\,,r 

= ||T(e-^g,)(T-i(l-XHH 



<^exp(-r^)||T-ni-XH) 



q,r 

v\ 



\q,r 



1 2n 

Using now the information that I T~ {x)\ < C e.xp{'y\x\^^) andthat(l —Xrv){x) = when 
\x\ > 4R/3, we get 



T-\l-XR)v\\^,r ^ C/exp(7(4i?/3)^)| 



so that finally 

Thus we have proven: 



< C exp (7(4i?/3)5^)exp(-T 



\q,r 1 



fc+ 1 

2n 



\t,m 



Proposition A. 4. Under the hypotheses of Proposition A.l, there exist constants C{(-^m) > 
and 7 > 0, such that for all v e H^ ,^, one has 



e^^'gf (1-Xfl)^ < Cexp(7(4P/3)^)exp(-r^)||^||,, 



k+1- 



We now return to the: 



Proof of Proposition 2.1. As in [EWW] it is only necessary to consider the term with 
highest derivative in 1 1 X^e^''^ 1 1 ^ • AH other terms are easier to estimate. Also, as in that paper, 
we use the fact that 

D^e^^ = e^^/^^^^e^^D^ , (A.7) 
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where is a shorthand notation for a product of derivatives w.r.t. the Xj of total degree £. Thus, 

rd/(2n) /• 

{e^^D'v){x) = ^-j^ J d''z9{z,T){D'v){e^/^'-\x + z)) . (A.8) 
First consider the q — case of (A. 3). Then 

Td/(2n) p 

IIXi^e-^^L,^ < y d''z\g{z,T)\\\w-^XRi^'v){e^/^'-\. + z))\\^, (A.9) 

where is the operator of muItipUcation by ( 1 + x • x) Note that the conclusions of Lemma 
B.4 of [EWW] do not depend on the exact form of g and so it also holds in the present situation 
and we have 



Lemma A. 5. One has the bounds 

\yXRv{e^^^'^\. + z))\\: < 



Ce"™ll^llo,m , if\z\<lR/S, 
2 - ^ C(l + \z\^y\\v\\l^, if\z\ > IR/S . 

Remark. Note that the proof in [EWW] is also unaffected by the dimension d in which we 
work. 

Now use Lemma A. 5 to bound the integral in (A.9) by writing it as an integral over 
\z\ <lR/% and an integral over \z\ > 7R/S. The integral over \z\ < 7i?/8 is bounded with the 
aid of Lemma A.5 as 

C 



grd/(2n) 



d 



(271) 



/ \9iz,r)\e--"'/'\\v\\,^^ < C{n)e'^'l^-e-^^\\v\\,^^ , (A.IO) 

J\z\<lR/8 



where the last step used the estimates of Theorem 4.1 to show that / dz\g{z, r) | < C, with C 
independent of r. 

To estimate the integral in the outer region, we use the second part of Lemma A.5 and then 
bound it by 



rd/(2n) /■ 



(A.U) 

for some 7 > 0, where, again, we have used the estimates of decay in Theorem 4.1 both 

— 7i?2n-i j as well as to bound the integral over z. Combining 
Eqs.(A.l), (A. 10), and (A. 11), we get the 9 = case of Eq.(A.3). 

We next indicate how to treat the g > cases of Eq.(A.3). Consider the case q = \. We 
can rewrite Eq.(A.8) by integrating by parts once w.r.t. one component of z, for example z-^. 
Then, 

rd/(2n) -r/(2n) . 

(e-^D^^)(a:) = j d''z{X^^^g){z,T){X^'-'v){e-/^^-\x + z)) . (A.12) 



Cahn- Milliard Equation 



16 



Differentiating (A.2) w.r.t. gives 

(D^^g){z,T) = J d'^kiki exp(z9 • z) exp((fc • - e"^)) . (A.13) 

To estimate (A.13), first replace k hy pj = a{Ty^^^'^^kj, where, as before a(T) = 1 — e~^. 
Then, 

i^.9)M = ^^^^;(,.) ^^^y,,,^, I d'^pPiexp((p.p)")exp(zp../a(r)V(^-)). (A.14) 

The estimate of the integral in Eq.(A. 14) now follows as before, since the extra factor of does 
not cause any trouble as it is easily offset by the exponentially decaying terms. 
One now uses the Schwarz inequality to rewrite 

rd/(2n) p 

WXRe'-^D'vlkm < J d'^z \D^M^,t)\ «;"^Xfi(D^-^) (e^/'(. + ^)) ^ , (A.15) 

and then proceeds as in the case when q = 0, breaking the integral over z into the same two 
pieces as before. These two pieces are then estimated with the aid of Lemma A. 5. Note that 
while the factor a(r)~'^/'^^"-* of Eq.(A.14) will be absorbed when one integrates w.r.t z, the 
remaining factor of a(r)-i/(2'') will remain in the final bound of Eq.(A.3). The bounds for 
Q = 2, 3, . . . 2n — 1 follow in a similar fashion. 

To complete the proof of Proposition 2. 1, first rewrite 

e-^Q, = e-^/^Q,e-^/^ = e^^/'Q.XRe^^/' + e^^/^Q,(l - X«)e-^/^ . (A.16) 

The second of these terms involves an estimate of the action of e^'~'^^Q^ on a function localized 
near the origin, so by Proposition A.4, we get a bound 

\V'-"Q.(^-Xny'-"<'h,m < C,exp(7(^)'"'<'"-" -i^j^) ll«ll,-„„ ■ (A.17) 
We use Proposition A.l to bound the first term of (A.16): 

^e^'/^(e-S™ + e-«'""-'>)||«||,_„„. (A.18) 
As a preliminary step, we note that if we first choose r and R such that 

grg/2 J^g-irr ^g-7fl^"/P"-')^ < ^- n{{k+\) / {2n)) ^ 

then for sufficiently small (roughly speaking ~ i (l + (4/3)^''/^^""^)) "\ the Eqs.(A.17) 
and (A.18) imply 

This shows that the projection of the semi-group onto the complement of the eigenspace 
spanned by the first k eigenvalues decays with a rate proportional to the eigenvalue A;,_|_j. We 
can sharpen the decay rate so that we obtain a rate like exp(— (1 — e) |A^^j |) by the techniques 
of [EWW], (see Eq. B.14 and following) and this completes the proof of Proposition 2.1 
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